function R=PlasticZone(KI)
global MatProp
Sy=365.422;
R=1/2/pi*(200/Sy)^2;